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Abstract 

A method is provided for approximating random slow manifolds of a class of slow- 
fast stochastic dynamical systems. Thus approximate, low dimensional, reduced slow 
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systems are obtained analytically in the case of sufficiently large slow-fast scale separa- 
tion. This dimension reduction procedure is applied to examine the impact of random 
environmental influences on the motion of inertial particles in a cellular flow field. The 
impact is partially quantified by computing first exit time and escape probability of 
inertial particles. 

Key Words: Random slow manifolds, dimension reduction, stochastic differential 
equations (SDEs), approximation under big scale-separation, inertial particles in flows 
Mathematics Subject Classifications (2010): 37H10, 37M99, 60H10 

1 Introduction 

Complex dynamical systems in science and engineering often involve multiple time scales, 
such as slow and fast time scales, as well as uncertainty caused by noisy fluctuations. For ex- 
ample, aerosol and pollutant particles, occur in various natural contexts (e.g., in atmosphere 
and ocean coasts [H HJ E|) and engineering systems (e.g. spray droplets), are described by 
coupled system of differential equations. Some particles move fast while others move slower, 
and they are usually subject to random influences, due to molecular diffusion, environmental 
fluctuations, or other small scale mechanisms that are not explicitly modeled [13]. Invariant 
manifolds are geometric structures in state space that help describe dynamical behaviors of 
dynamical systems. A slow manifold is a special invariant manifold, with an exponential 
attracting property and with the dimension the same as the number of slow variables. The 
reduced system on a slow manifold thus characterizes the long time dynamics in a lower 
dimensional setting, facilitating geometric and numerical investigation. 

Existence for slow manifolds of stochastic dynamical systems with slow-fast time scales 
has been investigated recently [T2l IT]. However, stochastic slow manifolds are difficult to de- 
pict or visualize. Therefore, in this paper, we approximate these random geometric invariant 
structures in the case of large time scale separation. We derive an asymptotic approxima- 
tion for these stochastic manifolds, and illustrate the random slow manifold reduction by 
considering the motion of aerosol particles in a random cellular fluid flow. The reduced slow 
system, being lower dimensional, facilitates our understanding of particle settling. 

We comment that approximations for individual solution paths (not a stochastic slow 
manifold) for stochastic slow-fast systems have been well investigated [El EJ [10] . Approxima- 
tions for deterministic slow manifolds have also been considered [9] [TTj. 

This paper is organized as follows. An approximation method for random slow manifolds 
is considered in £J2], and the dynamics of aerosol particles in a random flow field is investigated 
in P 
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2 Approximating random slow manifolds and dimen- 
sion reduction 



We first examine the existence of a random slow manifold for a slow-fast stochastic dynamical 
system, then devise an approximation method for this slow manifold, and thus obtain a low 
dimensional, reduced system for the evolution of slow dynamics. 

We consider the following slow-fast system of stochastic differential equations (SDEs) 

x = Ax + f{x,y), iGl", 
V = \By + lg(x, y) + fW t , y G P m . 

Here A and B are respectively n x n and mxm matrices. The functions / : M. n x M m — > M. n 
and g : P n x P m — > P m are C 1 -smooth and Lipschitz continuous with Lipschitz constants 
Lf and L g , respectively. The parameter a is a positive number and the parameter e > is 
small (representing scale separation). The stochastic process {W t : t G M} is a two-sided Un- 
valued Wiener process. When / and g are locally Lipschitz but the system has a bounded 
(i.e., in mean square norm) absorbing set or a random global attractor, a useful trick is to 
cut-off so that so that we still have a system with global Lipscitz nonlinearities and the new 
system has the same long time random dynamics as the original system. The existence of a 
random slow manifold for this system has been considered in [12] but we adopt a method in 
our earlier work [7]. 

We recall the definition of a random dynamical system (RDS) in a probability space 
(f2, J 7 , P). Let 9 = {8 t } t( z R be a B(K) <g) J 7 , J 7 — measureable flow, i.e., 

0:Rxfl4fl, 9o = id n , 9 tl o 9 t2 = 9 tl+t2 , fort^eR. 

Additionally, the measure P is supposed to be an invariant measure for 9 t , i.e. t F = P, for 
all £ G PL For a Wiener process driving system, we take Q = Co(R, R m ) consisting of all 
continuous sample paths of u(t) on P with values in P m and w(0) = 0. On Q, the flow 6 t is 
given by the Wiener shift 

9 t u(-) = w(- + 1) - w G f2, tel. 

A measurable map ip : M + x Q x P n+m — >■ M n+m is said to satisfy the cocycle property if 

<p(0, co, x) = x, (p(t + s, oj, x) = ip(t, 9 s oj, <p(s, uj, x)), for s, t G M + , oj G Q and x G P n+m . 

A random dynamical system consists of a driving system 9 and a measurable map with the 
cocycle property. 
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Introduce a Banach Space C\ as our working space for random slow manifolds. For A > 0, 
define 



C\ = \v : (—oo,0] — > R n : v is continuous and sup \e xt u(t)\ R n < oo}. 



and 



C\ — \y : (—oo,0] — > R m : v is continuous and sup \e xt v{€)\ Rm . < 00}, 

t<o 

with the following norms respectively 

\v(t)\ C { = sup \e xt v(t)\ R n, 

A i<0 

and 

|z/(t)| C 2 = sup \e xt u(t)\jim. 

Let C\ be the product Banach space C\ := C\ x C|, with the norm 

\(X,Y)\ Cx = \X\ ci + \Y\ cl . 

For matrices A and B, we make the following assumptions: 

HI: There are possive constants a, (3 and K, satisfmg —(3 < < a and K > 0, such 
that for every x G -R n and y € _R m , the following exponential estimates hold: 

\e At x\ R n < Ke at \x\Rn, t < 0; \e Bt y\ R m < Ke- pt \y\ Rm , t > 0. 
H2: > KL g . 

In order to use the random invariant manifold framework [5], we transfer a SDE system 
into a random differential equation (RDE) system. Introduce the following linear Langevin 
system 

dy = —ydt + -^=dW t . (2.2) 
It is known [5] that the following process ri e a (uj) is the stationary solution of the linear system 
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Moreover, 



Similarly, r] a (u) is the stationary solution of the following linear SDE system 

dy = Bydt + a dW t , (2.3) 

with 

/•O 

and 



7 ]a (0 t uj) = a j e Bit ~ s) dW s 4 <777(0t 



Denoting H^^a;) = 4^Wte(a/), it is also a Wiener Process [15] . and it has the same 
distribution as W t (u), with ijj e : — >■ Q, . Therefore, we have 

V £ (9 t£ uj) = ±= r e B(t " f } dW- = f e B ^ d{^W S£ {u)) = 7,(6^), (2.4) 



and 



if a {u) = 4= T e-^dW, = T e-^d(4^M) = JjftM- (2.5) 

In the calculation above, we used a transformation s' — s/e at the second equal sign and then 
omitted the prime in s'. Equations above indicate that r] £ (6tsUj) and r] £ (u) are identically 
distributed with rj(6 t i(} e oj) and 77 (^a;), respectively. 

We then introduce a random transformation 

X\ / x 

Y 



:=V g { U ,z,v)=l_ (nf{u) \, 



where (x, y) satisfies system (12.11) . 

Then the SDE system (12. ip is transferred into the following RDE system, 

X = AX + /(X, Y + arj £ (9 t u)), X e W 1 , 
Y = —Y + -g(X, Y + ar) e (6 t uj)), Y e R m . 



(2.6) 



By the variation of constants formula, this RDE system is further rewritten as 

X(t) = e tA X{0)+ f e A ^f(X(s),Y(s) + ar]%e s u;))ds, (2.7) 

Jo 

Y(t) = e B ^rY(t') + - [ e Bt -^g(X(s), Y(s) + a V £ (6 s u)) ds. (2.8) 

£ Jt' 

As y(-) G C\, we have the following estimation, 

\e Bt -^Y(t')\ Wn < e^\Y(t')\ Rm = e M ' \Y {t')^/ 1 ^ 1 0, f -oo, 
for A satisfying (3 — Xe > 0. 

Letting i! — > — oo, we get the expression of the RDE system (12. 6p . 

X{t) = e tA X(0)+ f e A ^f(X(s),Y(s) + ari £ (9 s uj))ds, (2.9) 

Jo 

Y{t) = - e Bt -^g{X{s),Y{s) + ar] £ {6 s u))ds. (2.10) 



We rescale the time by letting r — t/e, from system (12.61) and by (12.41) we get, 

X'(rs) = £[AX(re) + f(X(Ts),Y(rs) + ari(9 T jp £ u}))], X G M n , (2.11) 
Y'(re) = BY + g(X(re), Y{re) + ar](6 T ip £ uj)), Y G IR m , (2.12) 



where ' = 4-. 

dr 



We can rewrite these as the integral form below, 

X(re) = X(0)+e [ [AX(ss) + f(X(se),Y(se) + arj(9 s jp £ u;))] ds, (2.13) 

Jo 

Y{ts) = I e B{T - s) g(X(se),Y(se) + ar](6 s ij £ u))ds. (2.14) 

J — oo 

2.1 Dimension reduction via a random slow manifold 

We now recall some basic facts about random slow manifolds and dimension-reduced systems, 
when the scale separation is sufficiently large. 
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A random set M.(lo) = {(x, h(x,u))\x G IR n } is called a random slow manifold (a special 
random inertial manifold) for the system (12. ip . if it satisfies the following conditions [T2j : 

(i) Ai is invariant with respect to a random dynamical system 93, i.e. 

(p(t,u, M{u)) C M(6 t u) for t>0, ueQ. 

(ii) h(x,u) is globally Lipschitz in x for all uj G f2 and for any x G W 1 the mapping 
uj a;) is a random variable. 

(m) The distance of tp(t,u,z) and Ai(9 t uj) tends to with exponential rate, for z G R n+m , 
as t tends to infinite. 

A random slow manifold Ai, which is lower dimensional, retains the long time dynamics 
of the original system (12. ip . when e is sufficiently small [7]. 

In [12], a random Hadamard graph transform was used to estimate the existence of a 
random inertial manifolds, here we use Lyapunov- Perron method to achieve our result as in 

®- 

Lemma 1. Assume that HI and H2 hold and that there exists a A such that f3 — Xe > 0. 
Then, for sufficiently small e, there exists a random slow manifold Ai £ (u) = (£, h £ (^,u)) for 
the random slow-fast system ( 12. 6p . 

Proof. This proof is adapted from [7] for our finite dimensional setting. For completeness, 
we include the essential part here. For a A > 0, we use the Banach Space C\ as defined in 
the beginning of this section. 
Denote a nonlinear mapping 

T ((X, Y),X(0),u) = |e M X(0) + jT e A{t ^f(X(s), Y(s) + a^{6 s uj)) ds, 

- [ e Bt -^g(X(s), Y(s) + ar, e (6 s uj)) dsX. 
^ J —00 J 

Note that T is well-defined from R n+rn x W n x Q ->■ W n+rn . We will show that for every 
initial data X(0) = £ G M™, ( T276T) have a unique solution in C\. For (X, Y), (X, Y) G C A , we 
have that 
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T((X,Y),Z,u)-T((X,Y),Z,u> 

f e A ^ [f(X(s), Y(s) + af(9 s u)) - f(X(s),Y(s) + ar) £ {B s u))) ds 
Jo 



sup e 

t<o 



+ sup e 

t<o 



xt 



[g(X(s), Y(s) + arj%6 s uj)) - g(X(s), Y(s) + ar)%e s uj))] ds 



< 



(X, Y) - (X, Y) sup \e {x+a)t KL f [ e ~ {x+a)s ds) 
{X, Y) - (X, Y) sup { -e^-^KLg [ e ( l" A)s ds 
) (X,Y)-(X,Y) 



~~ ^a + A p-e\' 

The first inequality is by HI and the Lipschitz continuity of / and g, while the sec- 
ond inequality comes from direct calculation. Taking A = ^~^ Lg > 0, which satisfies 

f3 — Xe = P+^Lg > q ? we conc i uc l e that 



T((X 7 Y),Z,u))-T({X,Y),Z,u) 



< ( 



2KL f e 



2KL n 



2ea + (3-KL g KL q + f3' 



+ 1^T^) (X,Y)-(X,Y) 



By the assumption H2 



KL J1 y as e — > 0. Therefore, for £ small 



enou 



2KL f e 



2KL 



< 1, and 



+ XL qjg ) < 1. The contraction map theorem implies that for every 
£ G -R", T((X, Y), £, u;) has a fixed point (X(t), G C i3-kl 9 which is the unique solution 
of the differential equation system ( 12.6H . Moreover, the fixed point has the property 



(X(-;Z,u),Y(-,^u)) - (X(-;lu),Y(-,Z,u)) 
1 



< 



2KL f e 2KL g 
2ea+/3-KL g KL g +f. 



(2.15) 
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Denoting h £ (^,u) = Y(Q,£,u>), we obtain 

i r° 

h e (i,u>) = -l e- B ^g(X(s),Y(s) + ari £ (9 s u))ds, (Gl n (2.16) 
With the help of inequality (|2.15p . we further have 

2KL a K — £ ran 



h%Z,u)-h%£,uj) 



< 



KL g + 8 1 _ ( 2KL f £ + 2KL o ) 

y 1 \2ea+B-KL a ^ KL„+B ) 



Thus, h £ is Lipschitz continuous. By the fact that (X(0), Y(0)) G Ai e (u) if and only if there 
exists (X, Y) G C x and satisfies O and (12~TU|) . it follows that (X(0), Y(0)) G M £ (u) if and 
only if there exists such that (X(0),Y(0)) = (£, h E (£, oj)). Therefore, there exists a 

random slow manifold 

M°(u) = {{t,h*(M l^ffi"}. 

□ 

With the random slow manifold determined above, the dynamics on this manifold is 
described by the following dimension-reduced system in M. n (from equation (12. ip ), for e 
sufficiently small: 

£ = A£ + f(Z,<rv(M + h e (t,u)), eeK". (2.17) 
2.2 Approximation of a random slow manifold 

We now approximate the slow manifolds for sufficiently small e. Expand the solution of 
system (12. 12ft as 

Y(re) = Y (r) + eY x [r) + e 2 Y 2 (r) + ■ • • , (2.18) 
and the initial conditions as 

Y(0) = h £ (Z,u) = h^(^u) + eh^(^u) + ■■■, 
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and X(0) = (el n With the help of (T2~T3]) and ( l2~T8j) . we have the expansions 

f(X(TE), Y{re) + (rn{e T M) = /(£, Y (t) + ari^M) + /*(£, *o(r) + (t^tM) 

■e [ [AX(se) + f(X(se), Y(se) + a V {6 s ^ £ u))] ds 
Jo 

+fy(t,Y Q (T) +<tv(0tM)-[£Y 1 (t) + ■■■]+■■ ■ 

= m,Y o (T)+<TV(0rM) 

+m, Y {t) + ar){e T ^ E cu)) e [ [A£+ /(£, Y (s) + ^(0^))] da 

</o 

+/„(£, F (r) + t777(M^)) • eFi(r) + ■ ■ • , 
<7(X(re), K(re) + ar)(d T i/; £ uj)) = g(£, Y (r) + ar)(6 T i> £ u)) 

+£*(£, Yo(r) + vvWtM) -e [ [A£+ /(£, Y (s) + a V (0 s i> £ uj))] ds 

Jo 

+g v (£, Y (t) + arj(9 T jp £ oj)) ■ eY x {r) + ■ ■ ■ . 

Inserting ( 12. 18ft into ( 12.121) . expanding ( 12.12ft and then matching the terms of the same power 
of e, we get 

fy '(r) = BY (r) + g(C,Y (r)+^(e T M), (C) . 

\Yo{0) = W(Z, U ), ^ 

and 

f Y((r) = [B + g y {i,Y {r) + art^uj^Y^T) 

+g x (£, Y (t) + t777(M e w)){Are + fo Y (s) + a V (9 s ^ £ co)) ds}, (2.20) 

Solving the two equations for Yq(t) and Yi(r), we obtain 

Y {r)=e BT h (0 \i ) uj)+ [ T e- B ^g(Z,Y (s) + a V (9^ e uj)) ds, (2.21) 



o 
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and 

Y^T) = e BT+ fo By^YoisnvviOs^^ds^il)^^ ^ + f e ^B(s-r)+J^ g y (^Y (r)+a V (e r ^uj)) dr 

JO 

■g x (t,Y (s) + ar,(9 T i> e ou))[Ast+ [ (f(t,Y (r) + a V (9 r M)dr]ds. 

Jo 

(2.22) 

With the help of ( 12 . 1 3[) and (I2.18p . the expression (I2.16P can be calculated as follows 
h e {£,«>) = - f e- B ^g(X(s),Y(s) + ar ] £ (9 s u))ds 

^ J — oo 

r° 

/ e- Bs g(X(se), Y(se) + a V (8 s ij £ u)) ds 

J — oo 

g{£, Y Q (s) + a V (6 s iP £ u)) + g x (£, Y (s) + a V (e s ^ £ u))e [As£ 
+ I (/(£,y o (r) + a7K0^ 

e- Bs g(t,Y o (s) + ar)(0 s i/; e u;))ds 

oo 

+e f e- Bs \g x (t,Y {s) + <7V{8sM)[Ast+ [ (/(£, Y (r) + arj{e r ^ e u))) dr] 

J-oo JO 

+g y (£, Y (s) + ar ] (9 s ij £ u))Y 1 (s)} ds + 0(e 2 ). 



oo 




e~ Bs 



oo 

s 



e 2 ) 



To get the second equation, we used r = s/e and then used s to replace r. Thus the first 
and second order terms in e, of the random slow manifold Jvl £ (u) for ( 12.61) . are respectively 

h^^,u)= f e- Bs g^,Y o (s) + ar](0^ e cj))ds, (2.23) 

J — oo 



and 



hW(£,u) = f e~ Bs \g x (t,Yo(s) + ar)(e s i> e cj))[Ast+ [ /(£, Y (r) + a V (0 r ^ e cu)) dr] 

J -oo L JO 

+^(e,io(s) + <r»7(0 s ^))yi(s)}ds. (2.24) 
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That is, the slow manifolds A4 £ (u) = {(£, h £ (!;,u))} of (12. 6p up to the order 0(e 2 ) is repre- 
sented with h £ (^,u) = h(°'(£,oj) + eh^(^,co) + 0(e 2 ). This produces an approximation of 
the random slow manifold. 

Therefore, we have the following result. 

Theorem 1 (Approximation of a random slow manifold). Assume that HI and H2 hold and 
assume that there is a X such that (3—\e > 0. Then, for sufficiently small e, there exists a slow 
manifold M £ ((d) = {(£,h e (£,u) | £ G M n }, where h £ (^ u) = h (0) (£,a/) + eh w (^id) + 0(e 2 ) 
with h(°\t; , u) and h^ 1 '^, u) expressed in (12. 23ft and (12. 24ft . respectively. 

With the approximated random slow manifold 

h £ (^u)= h^(^oj) + eh^^u), (2.25) 

we obtain the following dimension-reduced approximate random system in MJ 1 (from equation 
(I2.17P ), for e sufficiently small: 

£ = AZ + f^ar ] (ij £ id) + h £ ^id)), eeR". (2.26) 



3 Settling of inertial particles under random influences 

For the motion of aerosol particles in a cellular flow field, Stommel once observed that, 
ignoring particle inertial (e = 0), some particles follow closed paths and are permanently 
suspended in the flow. Rubin, Jones and Maxey [11] showed that any small amount inertial 
(small e > 0) will cause almost all particles to settle. Via a singular perturbation theory [9], 
Jones showed existence of an attracting slow manifold. By analyzing the equations of motion 
on the slow manifold, especially heteroclinic orbits, they could understand some mechanisms 
that inhibit trapping and enhance settling of particles. 

Let us now examine the motion of aerosol particles in a random cellular flow field, using 
the slow manifold reduction technique developed in the previous section. 

Consider a model for the motion of aerosol particles in a cellular flow field, under random 
environmental influences [TTJ 



yi = vi, 
m = v 2 , 

vi = + \a sin?/! cos y 2 + ^W?, 
V2 = -l V2 + ±{y - a cos yi sin y 2 ) + 
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where (2/1,2/2) and (i>i, 1*2) are position and velocity, respectively, of a particle in the horizontal- 
vertical plane (positive 2/2 axis points to the settling/gravitational direction), a is a velocity 
scale, and V is the settling velocity in still fluid. Moreover, W t , Wf are independent scalar 
Wiener process, a is a positive parameter, and e is the inertial response time scale of the 
particle. Note that (a sin 2/1 cos 2/2, —a cos 2/1 sin 2/2) is the so-called cellular flow field on the 
domain (a 'cell') D = (0, vr) x (0,vr). 

As in Section [2j this four dimensional SDE system can be converted to the following 
RDE system 



Denoting 2/i(0) = £1 and 2/2(0) = £2, and we examine the motion of the particle (£1,^2)- By 
using (123511 and (j272T|) . we get 



'ill 



v 1 + ar]f(9 t uj), 

-\vi + ~a sin 2/1 cos 2/2, 
--^2 + ~(V - a cos 2/1 sin 2/2), 




(3.2) 



U'2 



where 








e s a sin ^ cos ^ 2 ds = a sin £ x cos £ 2 , 



/4°^ (£, w) = V — a cos £1 sin £ 



—00 



and 



Uio(£) 



asin£i cos£ 



2, ^2o(^) = V — acos£isin£ 2 - 
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Owing to f[2~24l) . 

e s [a 2 s sin£i cos£x — aVs sin£i sin£ 2 + a cos£i cos £ 2 a / r]i(6 r t/j £ uj) dr 

-oo Jo 

— a sin^ sin ^ 2 0" / V2(0 r i/j £ u) dr\ ds 
Jo 



—a sin^i cos^i + aV sin£i sin£ 2 + aer cos £i cos £2 / se s dW s 
— a<7sin£isin£ 2 / se s dW 2 





00 



^2^(^15 ^2) ^) — —a 2 sin ^2 cos ^ 2 + cos ^1 cos ^ 2 + sin ^ sin ^2 / se s dW] 

J — 00 

/•O 

— aa cos £1 cos £2 / se s dW 2 . 



Therefore, from (12.171) . the dynamics on the random slow manifold is described by the 
following dimension-reduced system: 

6 = ^1(61, £2, w) + <nji tysu) 
f° 

= a e s dW s (il> £ u) + a sin£i cos£2 + s{ — o 2 sin£i cos£i + aV sin£i sin£ 2 
J —00 

/0 pO 
se s dW] — aa sin£i sin£ 2 / se s dW 2 }, (3.3) 
-00 J — 00 

6 = M£i,6,w) + arj 2 (i> £ u) 

— a I e s dWg(ip e uj) + — a cos£i sin£ 2 + ^{ — « 2 sin£ 2 cos £2 + aV cos£i cos £2 
J —00 

/0 /-O 
se s dW] - aacos^i cos£ 2 / se s dW 2 }. (3.4) 
-00 j —00 
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3.1 Numerical simulation: First exit time and escape probability 

Note that (£1, £2) is the particle position. For random slow manifold reduction, it is customary 
to use a notation different from the original one (2/1, 3/2)- The positive £2 direction points to 
the bottom of the fluid. 

In this section, we conduct numerical simulations for this reduced or slow system (13.31) - 
(13 .4 j) . When e = 0, a = 0, this reduced system becomes the system for the motion of 
particles in the cellular flow. When e 7^ 0, a = indicates no noise, while a non-zero a 
means noise is present. 

Simulation Motivated by understanding the settling of particles as in [11] . we first cal- 
culate first exit time of particles, described by the random system ( I3.3p -( l3~4l) . from domain 
D — (0, 7r) x (0, 7r) and then examine how particles, exit or escape from the fluid domain D. 
To this end, we introduce two concepts: First exit time and escape probability. The first 
exit time is the time when a particle, initially at (£1,^2) £ D, first exits the domain D: 

r(£i,&)=krf{*: (£i(t)Mt)) t D}. 

Let T C dD be a subboundary. The escape probability P r (£ 1 ,£ 2 ), for a particle initially at 
(£1, £2), through a subboundary T, is the likelihood that this particle first escapes the domain 
D by passing through T. We will take T to be one of the four sides of the fluid domain D. 
The escape probability of a particle through the top side V = {£ 2 = means the likelihood 
that this particle settles to the bottom of the fluid (note that the positive £ 2 direction points 
to the bottom of the fluid). 

To compute the first exit time from the domain D, we place particles on a lattice of 
grid points in D and on its boundary, and set a big enough threshold time T. As soon as a 
particle reaches boundary of D, it is regarded as 'having exited' from D. If a particle leaves 
D before T, then the time of leaving is taken as the first exit time, but if it is still in the 
domain at time T, we take T as the first exit time. When a particle's first exit time is T, we 
can see it as trapped in the cell. 

In order to calculate the escape probability of a particle under noise through a subbound- 
ary T, one of the four sides of the domain, we calculate a large number, N, of paths for each 
particle to see how many (say M) of them exit through the subboundary T, and then we 
get the escape probability We do this for particles placed on a lattice of grid points in 
D and on its boundary. When a particle reaches or is on a side subboundary, it is regarded 
as 'having escaped through' that part of the boundary. 
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Figure 1: Phase portraits of the slow system (I3.3l) -( l3^4l) with V = 0, a = 0.7 and a = (no 
noise): e = (top); e = 0.05 (bottom). 
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Figure 2: Sample particle orbits of the slow system f l3 . 3 H and f 13.4f) with = 0.3 and a = 0.7: 
£ = and a = (top, no noise); £ = 0.05 and a = (middle, no noise); £ = 0.05 and a = 0.01 
(bottom, with noise). 
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Figure 4: Escape probability for the slow system (13.31) and (13. 4p with V = 0, e = 0.05, 
a = 0.7 and o = 0.01: (a) escape through £ 2 = tt (settling direction or physical bottom 
boundary), (b) escape through £2 = (physical top boundary), (c) escape through £1 = 
(left boundary), and (d) escape through £1 = it (right boundary). 
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Figure 5: Escape probability for the slow system (\3.3\i and H3.4[) with V — 0.1, e — 0.05, 
a = 0.7 and a = 0.01: (a) Escape through £ 2 = 7r (settling direction or physical bottom 
boundary), (b) escape through £2 = (physical top boundary), (c) escape through £1 = 
(left boundary), and (d) escape through £1 = it (right boundary). 
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Figure 6: Escape probability for the slow system ( 13. 31) and ( 13. 41) with = 0.1, e — 0.05, 
a = 0.7 and a = 0.01: (al) and (a2) are Figure[5](a) splitting at the deterministic heteroclinic 
orbit asin£i sin £2 — V£\ = 0; (cl) and (c2) are Figure M^c) splitting at the same heteroclinic 
orbit. 
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Figure 7: Deterministic stable manifold H^ s (0, 7r — sin 1 (^-)) (Blue solid curve) and unstable 
manifold W u (0, sin _1 (^)) (Blue dashed curve) for the slow system ( 13.31) and (13 .4p : V = 0.1, 
e = 0.05, a = 0.7, and a = (no noise). Two equilibrium points are (0,7r — sin _1 (— )) and 

(O.sin" 1 ©)- 

3.2 Deterministic Case 

We recall the results, by one of the present authors, of the deterministic motion of aerosol 
particles in a cellular flow [TT] . 

When the settling velocity in still fluid is zero, i.e., V = (and also a — 0), particles are 
trapped in the cell either in circular motion (with no inertial, e = 0) or spiralling motion 
(with inertial, e > 0), as shown in Figure [T] (top) and (middle), respectively. 

When the settling velocity in still fluid is non-zero, i.e., V > 0, in the case with inertial 
absent (e = 0) and noise absent (a = 0), the particles in the area surrounded by the 
heteroclinic orbit a sin ^ sin £ 2 — V£i = connecting the equilibrium points (0, sin _1 (V/a)) 
and (0, 7T — sin^ 1 (V / a)) , are trapped inside it, with the equlibrium point (cos~ 1 (K/a), n/2) as 
a center. But the particles in the remaining area settle to the bottom of the fluid. With an 
arbitrarily small inertial effect (0 < e 1 and also a = 0), the heteroclinic orbit breaks and 
it leads to the settling of almost all particles, with the equilibrium point (cos _1 (\^/a)), 7r/2) 
becoming an unstable spiral. Figure [7] is the stable manifold ^(0,^ — sin _1 (\^/a)) (Blue 
solid curve) and unstable manifold If" (0, sin _1 (y/a)) (Blue dashed curve) when inertial 
presents (e = 0.05 and also a — 0). 
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Figure 8: Escape probability (Light blue color) for particles through the left boundary 
for the slow system (E3D and (E3D with V = 0.1, e = 0.05, a = 0.7 and a = 0.01 (with 
noise), together with deterministic stable manifold W^ s (0, 7r— sin _1 (— )) and unstable manifold 
W u (0,sin _1 (^)) (Red curves). 
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3.3 Stochastic Case 



For zero settling velocity in still fluid (V = 0), we note that all particles are trapped in a 
fluid cell D = (0,7r) x (0, ir) when noise is absent. Figure [1] shows the particle orbits in this 
case. 

For non-zero settling velocity in still fluid (V ^ 0), when noise is absent (er = ), all 
particles settle to the fluid bottom; see Figure [2] (top, middle). But when noise is present 
(er 7^ 0), some particles exit the cell. Figure |2] (bottom) shows that, with a small noise 
(er = 0.01), some particles indeed exit the cell from the vertical side boundary £i = 0. 

In fact, when noise is present, all particles will exit, no matter the settling velocity in still 
fluid V is zero or non-zero. Figure [3] indicates that with noise, particles will all exit from a 
fluid cell in finite time. In the following we only consider the case with noise. 

Figures 0] — |5] plot the escape probability through four side boundaries, for zero or 
non-zero V (settling velocity in still fluid) values. When a particle reaches or is on a side 
boundary, it is regarded as 'having escaped through' that part of the boundary. In other 
words, particles on a side boundary have escape probability 1 (you see this in these figures). 
When V — 0, the particles escape the cell through each of the four side boundaries with 
similar or equal likelihood (Figure HJ), as there is no preferred direction for particles due to 
zero still fluid velocity V = 0. With non-zero V, particles almost not escape through the 
right side boundary £i = n. In fact, the inertial particles either settle to the physical bottom 
or exit from the left side boundary £i = 0. Figure [5] displays the escape probability for 
V = 0.1, through each of the four side boundaries of the fluid cell. Although most particles 
settle (Figure [5] (a)), some particles escape the fluid cell through the left side boundary 
(Figure El (c)). See Figure |6] for a split view of this phenomenon. 

To examine this phenomenon more carefully, we draw the stable manifold W s (0,n — 
sin _1 (V A /a)) and unstable manifold W u (0,sm~ 1 (V/a)) for the deterministic system (<r = 0) 
in Figure [7J As shown in Figure El inertial particles with significant likelihood escaping 
through the left side boundary are near or on the stable manifold VF s (0, ir — sin _1 (y/a)). In 
other words, inertial particles near or on this stable manifold are resistent to settling in the 
stochastic case. This resistance is quantified by the escape probability for a particle to get 
out of the fluid cell through the left side boundary. 

3.4 Conclusions 

Let the settling velocity in still fluid be non-zero (i.e., V ^ 0). 

(i) In the classical case (no inertial: e = and no noise: a = 0), the particles surrounded 
inside a heteroclinic orbit are trapped inside it and all the other particles settle to the bottom 



24 



of this cellular fluid flow. 

(ii) In the case with only small inertia influence (0 < e <C 1, o = 0), the heteroclinic orbit 
breaks up to form a stable manifold W s and an unstable manifold W u , the trapped particles 
then settle, i.e., all inertial particles settle. 

(hi) However, when the noise is present (0 < e -C 1, o~ ^ 0), although most inertial particles 
still settle, particles near or on the deterministic stable manifold W s escape the fluid cell 
through the left side boundary, with non-negligible likelihood. The deterministic stable 
manifold W s acts as an agent against settling. 
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